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We consider one particle confined to a deformed one-dimensional wire. The quantum mechanical 
equivalent of the classical problem is not uniquely defined. We describe several possible hamiltonians 
and corresponding solutions for a finite wire with fixed endpoints and non-vanishing curvature. We 
compute and compare the disparate eigenvalues and eigenfunctions obtained from different quanti¬ 
zation prescriptions. The JWKB approximation without potential leads precisely to the square well 
spectrum and the coordinate dependent stretched or compressed box related eigenfunctions. The 
geometric potential arising from an adiabatic expansion in terms of curvature is at best only valid 
for very small curvature. 

PACS numbers: 03.65.Ca,03.65.Ge,67.85.-d 


I. INTRODUCTION 

Advances in optical trapping techniques makes it pos¬ 
sible to trap cold atoms in what is effectively 1 or 2 
dimensional [il-Q. An interesting direction of this is the 
ability to create effective quantum wires of varying shape 
and curvature. An advantage of these lower dimensional 
setups are the how they can be shaped in ways that 
provide stabillity against short range collapse of polar 
molecules. Another thing is how long range interactions 
can shortcut through the second or third dimensioi^ and 
create novel interactions in 1 and 2 dimension [8l4l2|. 

When describing the motion of an atom moving along 
such a wire there are at least two different approaches. 
The first is to solve the problem fully in three dimensions, 
and treat the trapping potential as any other potential. 
This approach is difficult to solve both analytically and 
numerically, unless the trapping potential is highly sym¬ 
metric. The other approach is to develop a quantum 
mechanical description of the system in one dimension. 
How to do this is not well defined for general structures. 
This paper will discuss various quantization prescriptions 
of the classical motion, both in general and in particular 
focussed on the one-dimensional wire problem. 

In general, position dependent masses present prob¬ 
lems in quantizing classical motion. This problem was al¬ 
ready recognized by Schrodinger [T^ , and in the dynami¬ 
cal evolution of nuclear shapes from equilibrium through 
saddle points to separated fission fragments M- The 
problem appears first of all when the starting points are 
generalized masses obtained independent from potential 
energies as for example through classical cranking type 
of models [l5l| and micro-macro energy calculations [iq 

The motion can be constrained geometrically by an ex¬ 
ternal field and it is then possible to simulate this physical 
situation by steep confining walls in forbidden directions. 
Then the procedure can be more straightforward by ap¬ 
proximately solving the full many-dimensional problem. 
The effects of the repulsive walls may then be approxi¬ 
mated by applying an adiabatic expansion which results 
in a potential reflecting the geometric properties of the 
confinement 0. We shall compare results from different 


quantization prescriptions. 

The purpose of the present paper is to exhibit similari¬ 
ties and differences between quantum mechanical results 
arising from different quantizations of a single-particle 
hamiltonian in one dimension. In general we shall allow 
an undetermined coordinate dependent parametrization 
of the particle mass. In section m we develop a classi¬ 
cal description of a particle confined to move in along a 
curved wire. In section uni we show the different choices 
of how to quantize the system. In the sections IIVI and m 
we look at three different choices of wire, and using the 
different approaches from section Ell we calculate the 
spectra and the eigenfunctions. Finally in section |VT] we 
give a summary of the discussions and an outlook. 

II. CLASSICAL DESCRIPTION 

We consider one particle confined to move in a one¬ 
dimensional trap. The particle is assumed to be point 
like with no internal structure and have a mass of rriQ. 
The trapping geometry is going to be several different 
deformations of a helix. This one dimensional curve is 
defined by a parametrization A : K. —>■ R^. The parame- 
terizations discussed in this paper will all be of the form 

{x,y,z} = R-. (1) 

Instead of using the full three dimensional set of coordi¬ 
nates we use (j), to describe the position of the particle 
along the curve. 

The classical velocity is then {x,y,z), where the dots 
denote time derivation. The corresponding classical ki¬ 
netic energy, T, of the particle in three dimensions is then 
given by 

T = imo + i^) ■ (2) 

Then T can be transformed through Eq. ([T]) to only de¬ 
pend on the angular position, (j), and velocity, (j), along 
the curve, that is 

<P) = ^moix'"^ + y''^ + , (3) 
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which defines the effective mass For the 

parametrization in Eq. m this yields the explicit expres¬ 
sion 

m{<j3) = moR^ [ifj + /2) cos^ 4>+ifx+ f'y) sin^ 

+ (/./^-/y/y)sin(2</.) + /'2] , 

where the primes denote derivatives with respect to (f). 

A useful quantity is the canonical conjugate momen¬ 
tum, to the position cj) along the curve, that is [l^ 

P4> = ^= , (5) 

ocp 


operator and the mass term, l/m((^), is now important. 
We demand that the hamiltonian is hermitian and the 
eigenenergies of the system must consequently be real. 

A general structure of the kinetic energy operator is 

(9) 

2 m“ d(j) mP d(j) ’ 

where a -|- 6 -I- c = 1 and hermiticity requires a = c. An 
even more general form would be to allow a ^ but 
instead symmetrizing afterwards, that is 

P 1 d 1 d 1 1919l\ 

4 \ m“ d(j) mP d(f> mP d4> tvP dcj) m°-) 


which provides the kinetic energy in canonical form 


pI ^ 

2m{4>) 


( 6 ) 


The coordinate dependent mass in Eq. o is a key quan¬ 
tity. The present parametrization in Eq. © only pro¬ 
vides ^-independent effective mass when the functions 
/' and fx = fy are independent of (p. These are both 
necessary and sufficient conditions for the constant effec¬ 
tive mass, m{cj)) = -|- /'^), of a regular helix. 

Another simple case is fx = fy where we have 


m(<^)=moi?2(/J + /^ + /f), (7) 


which still may contain a (/)-dependence. When both 
fx and fy are constants, which corresponds to the 
parametrization of an ellipse in the x-y plane, we have 


The quantization in Eq. (113 is problematic for discon¬ 
tinuous potentials, which breaks the hermitian property 
M- This does not apply to our case and we shall not 
a priori reject the two different terms in Eq. (flUl) . To 
be specific we choose typical examples with properties 
from Eqs. ([9]) and (ITO)) . that is a = c = 0, &=1 and 
6 = c = 0,a = l, respectively. 

The first of these choices of an effective mass hamilto¬ 
nian is with the mass placed between the two momentum 
operators, that is 


d 1 d 

2 d(j) rn{cj)) dp 


( 11 ) 


By taking the first order derivative of the effective mass 
the hamiltonian can be rewritten as 


3^ h^m'{p) d 

2m{p) dp"^ 2m'^{p) dp 


( 12 ) 


m{p) = moR^ {fl cos^ P^ fl sin^ P + f?) ■ (8) 

If furthermore, fx = fy and /' = 0 the helix reduces to 
a circle in the x — y plane, and the mass becomes moR^, 
which is just the initial mass scaled by the i?^-factor that 
appears when using the dimensionless (/)-coordinate. We 
conclude that the effects of a non-constant effective mass 
require more than the regular helix, and we need to con¬ 
sider angle-dependent /-functions. 


III. QUANTIZING THE MOTION 

In this section we shall consider different prescriptions 
to quantize the particle motion on the curves parameter¬ 
ized through the classical physics described in the previ¬ 
ous section. 


A. Quantizing in one dimension 

Without an external potential on the particle the 
hamiltonian only contains terms arising from the kinetic 
energy in Eq. ([5]). The complication may be the position 
dependent effective mass. The straightforward quantiza¬ 
tion is from Eq. ([5]) with but ordering of this 


We note here that the position dependent effective mass, 
results in an extra term in the hamiltonian. Not only 
does this hamiltonian contain a first order derivative of 
the effective mass, but this term also contains a first order 
derivative. Both features are not present in a normal 
hamiltonian only containing a kinetic energy part with a 
constant mass. 

The first order derivative can be removed as usual by 
the substitution of the total wavefunction, T = 
into the Schrbdinger equation, (Hi — E)'^ = 0. The cor¬ 
responding new hamiltonian acting on the reduced wave- 
function, p then becomes 

2m{p) J ’ 
(13) 

where the energy remains unchanged while a centrifu¬ 
gal barrier term appears, and the Schrbdinger equation 
becomes {Hemi — E)p = 0. 

Another possible choice of hamiltonian is one explicitly 
constructed as hermitian with two terms. One term with 
the effective mass before the momentum operators, and 
another term with the effective mass after the momentum 
operators, that is 

tP f d^ 1 1 d^ 

^ 4 'm{p) m{p) dp'^ 



HeMi = — 


d^ 


2m{p) dp"^ 


4m^((/)) \ 


m"{p) 
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Calculating the derivatives of we end up with a 

hamiltonian, with the same first order derivative term 
appearing in Hi, and in addition another effective po¬ 
tential. In total we get 

tirin'{jp) d / 1 \ 

^ 2m{(j)) d(j)^ 2m?{(j)) dp 4 dp'^ \m{p) J 

(15) 

The extra potential term is the only difference between 
Hi and H 2 . It contains terms of up to second order in the 
derivative of the mass. The first order angular derivative 
is removed as in Hi resulting in the hamiltonian, Hem 2 i 
acting on the reduced wavefunction, p. The difference 
between these two hamiltonians amounts to an additional 
potential, that is 


Vem2 


HeMi — HeM2 

d^ r 1 \ 

4 902 \to(0) / 


(16) 



Thus, if the inverse of the effective mass by chance or by 
choice has a vanishing second derivative, the potential, 
Vem 2 : in Eq. (fT6|) is zero, and the two different schemes 
for quantization yields exactly the same results. 

Specifically, this means that Hemi = He M 2 if n given 
parametrization produces an effective mass such that 
l/m(0) = Co -f cip, where cq and ci are constants. We 
shall later study a system with such an effective mass. 
Extending to a quadratic power dependence in l/m{p) 
would produce a constant difference between Hemi and 
Hem 2 without leading to any difference in structure. We 
emphasize that the two choices of hamiltonian clearly 
only are examples. 


B. Quantizing in three dimensions 


A different approach to a quantum mechanical descrip¬ 
tion of a curved wire is found in [i3- They start with a 
potential in three dimensions, assume a very tight con¬ 
finement in the directions perpendicular to the wire, and 
apply an adiabatic expansion. Only the lowest excited 
transverse state is allowed and a one dimensional hamil¬ 
tonian is then derived, that is 


tP 92 ^2^2 

2mo 9s2 Sttio ’ 


(17) 


where s is the arc length along the wire, and k is the cur¬ 
vature of the wire as defined in the Frenet-Serret appa¬ 
ratus [13 . The last term in Eq. m is the geometric po¬ 
tential, and the curvature for a curve in parametrized 
by Eq. © is explicitly given in [13 to be 


2 {z”y' - y’'z'f + {x"z' - z"x'f + (y"x' - x"y'f 


(a;'2 -I- y'2 -I- 


(18) 

where the primes again denote derivatives with respect 
to p. 


The geometric potential is always attractive, and most 
attractive where the curvature is large. The curve is here 
parameterized by the arc length s, whereas we above used 
the azimuthal angle, p, to specify the position. The con¬ 
nection can be found by calculating s as function of p, 
and if necessary invert the resulting expression to get 
0(s). 

The arc length of a curve in defined by Eq. o can 
be calculated by [13 



(19) 


where we measured from p = pmin and used Eq. © 
and the related coordin ate depen dence on the angle, p. 
This means that s' = i/m((p)Jrno, which can be used to 
transform the hamiltonian in Eq. (ED from the s to the 
p coordinate. The result is 

hP 9^ fPm'{p) 9 
® 2m{p) dp"^ 2nn?{p) dp 

h?m"{p) 7h?m''^{p) tPtP 

8m2(0) 32m^{p) 8mo 

which can be verified most easily by going from Eq. o 
to Eq. (ITTI) by use of d/dp = ds/dp x d/ds. 

The geometric potential or rather the curvature, n, can 
be rewritten in terms of to, m' and to", and some extra 
terms containing third derivatives with respect to p, that 
is 


2 iTOTO" - iTO'2 - m{x'x'" + y'y'" + z'z'") ^ 

At — - o - , (21) 

TO'^ 

which can be verified by use of Eqs. © and (IT^ . Thus 
the curvature can be expressed through to, to', to", and 
additional terms containing more than third derivatives 
with respect to p. Again we find the same first order 
derivative term as in Hi and 772 • Removal results in the 
hamiltonian, Hgeo, acting on the reduced wavefunction. 

We now have three different hamiltonians describ¬ 
ing the same system. They are all one-dimensional 
Schrodinger equations with first and second derivatives 
as well as various potential terms. The first two, Hemi 
and Hem 2 ^ only have kinetic energy terms whereas the 
last one, Hgeo, also contains the attractive geometric po¬ 
tential. However, depending on quantization prescrip¬ 
tion, the kinetic energy parts differ from each other in 
these three cases, although all of them have the ordinary 
second order derivative term with the same coordinate 
dependent mass, that is — The same first order 

derivative term, ^, has been removed in all three 

cases. 

The hamiltonian, Hgeo, with the geometric potential 
differs substantially from the other ones. We note that 
the first two terms in Eq. ED containing derivatives of 
the mass are also present in He M 2 , except that they ap¬ 
pear with different strengths. However, the last term 





















4 


in Eq. (1211) has a different structure with higher order 
derivatives of the coordinates in three dimensions. 

We emphasize that the geometric hamiltonian is de¬ 
rived under various assumptions, and especially the adi¬ 
abatic approximation which cannot accommodate too 
rapid changes of the coordinates. This means in par¬ 
ticular that at most terms up to second order deriva¬ 
tives are correctly included, while derivative terms of 
third an higher order have been neglected. The accu¬ 
racy of such approximations is dubious, because higher 
order derivatives of the periodic trigonometric functions 
in the parametrization are not decreasing in size. This 
does not prove that the result is inaccurate but more in¬ 
formation is necessary to evaluate the consequences of 
these assumptions. 

The curvature is an important quantity, at least as 
long as it remains modest in size. It is then of interest to 
know when it vanishes, or equivalently when it becomes 
small. From Eq. (1181) we see that k = 0 is obtained when 
z"y' - y"z' = x"y' - x'y" = z"x' - x"z' = 0. These dif¬ 
ferential equations can be rewritten (Incc')' = (Inj/')' = 
(Inz')' with the complete solutions z = ay+ b = cx + d 
for arbitrary constants a, b, c, d. Thus, a linear depen¬ 
dence between all the coordinates eliminates the curva¬ 
ture. This is of course not surprising since a straight line 
by definition should have curvature zero. Still, it empha¬ 
sizes the point that only a modest correlated coordinate 
variation is allowed to maintain a small curvature and a 
fairly accurate adiabatic expansion. A helix seems to be 
far away from this assumption. 


C. Semi-classical approach 


The search for an appropriate quantization prescrip¬ 
tion of a classically well defined problem strongly sug¬ 
gest use of semi-classical methods. We therefore turn to 
the JWKB (Jeffreys-Wentzel-Kramers-Brillouin) approx¬ 
imation which directly is applicable on a one-dimensional 
problem. The lowest order expression for the wave func¬ 
tion of a bound state is 


'^jwkb{(I)) = A cos 


-I- B sin 


f {E-V{(p))d(p 

4*min 

( 22 ) 


where E is the energy of the particle, A and B are con¬ 
stants, 4’min is one end of the wire, and V{(j)) is the po¬ 
tential along the wire. The integral has to extend over all 
classically allowed regions of (p, that is where E > V{(j)). 
The expression in the exponent is in fact found as an inte¬ 
gral over the classical canonical momentum, derived 
from Eq. and the constraint of energy conservation 
T = E-V. 

This choice is not unambiguous due to the coordinate 
dependence of the effective mass. We could choose one 


of the hamiltonians as the starting point, then rewrite 
as in Eq. m where the first order derivative is removed 
and a reduced equation obtained. This is similar to the 
use of spherical coordinates and the equation for the re¬ 
duced radial wave function. Then the extra centrifugal 
potential should be included. We could also start with 
He M 2 s-nd only include the potential in Eq. (HU), or for 
that matter any linear combination of these potentials. 

The simple JWKB wave function in Eq. (|2^ could 
be improved but the full solutions are easily available 
to us for the different hamiltonians. We shall therefore 
only use the JWKB approximation to gain qualitative 
insight. It is obvious that both JWKB approximation 
and the geometric potential are only reliable for modestly 
varying coordinates along the wire and in turn slowly 
varying effective mass. The different hamiltonians are 
then also rather similar as their differences stem from 
the derivatives of mass or coordinates. Thus, it suffice to 
use E = 0 in Eq. (1^ . 

The boundary conditions we choose are precisely van¬ 
ishing wave function at the points terminating the classi¬ 
cally allowed regions as for example at the two ends of the 
finite wire. This immediately requires that A = 0 since 
'^.]WKb{4> = 4>min) = 0. The other end point condition 
of ^jwkb{4> = 4>max) = 0 then provides the general 
quantization condition, that is 


m{(j))d(j) = mr , (23) 

This equation is only fulfilled for discrete values oi E = 
En, and the corresponding wavefunctions are then given 
by Eq. (E^ for A = 0. 

The integral over is measuring the total length 

of the curve as seen by Eq mi). The spectrum is therefore 
exactly that of a particle in an infinitely deep square well, 
that is 

= d- T2J2 -V’ Lip) = / ^/m{(j))/mod(j) ,(24) 

2moL^[(j)rnax) J 

where L{<pmax) is the length of the one-dimensional box. 
The corresponding eigenfunction is 

'^JWKB oc sin (mrL{(j))/{L{(j) max )) • (25) 

We can now use this simple expression to compute the 
JWKB wave function for particular parameterizations 
where can be analytically integrated. Thus, the 

spectrum is that of a square well, and the related eigen¬ 
functions are deformed (stretched or contracted) one¬ 
dimensional box wavefunctions. 



IV. BULGING HELIX 

We shall now calculate the properties of the quantized 
structures. We must then first decide on an appropriate 
parametrization of the one-dimensional curve. Second 
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FIG. 1: A schematic drawing of a helix with a gaussian defor¬ 
mation of the type in Ea. (l26ll . The parameters are chosen to 
be a = 0.5 and and the radius then varies by a factor 

of 1.5 over about 4 windings. The length can be extended as 
desired. 

we compare numerically the resulting coordinate depen¬ 
dent masses, curvatures, spectra and eigenfunctions for 
different quantization prescriptions. 

A. Parametrization 

In section El we showed that a regular helix leads to a 
constant effective mass. The curvature is then naturally 
also small or at least varying slowly with position along 
the one-dimensional curve. The different quantization 
prescriptions are then identical, except perhaps for the 
one-dimensional adiabatic approximation of the three- 
dimensional result. To emphasize the effects of varying 
mass and curvature we therefore start by modifying the 
simple helix. We keep the circular nature, that is the ra¬ 
dial variation of the x and y-directions are chosen to be 
identical, but varying quickly over a few windings. The 
unrelated z-direction is linear to maintain the equidistant 
helix structure. In total, we choose the parametrization 
in Eq. © to be 

fxW = = (1 + a • exp(-(^V(^2)^ ^ 

= (t> , 

where a and tpo are constants. For practical convenience 
we use the gaussian to describe the form of the variation. 
The size of a determines the radial changes, ranging from 
R to i?(l-|-a) and back again as 4> varies from —oo to -boo. 
The width of the gaussian, (/)o, determines how quickly 
the radial change is taking place, that is over how many 
windings. The length of the curve is chosen to be from 
4>min = —40 to <t>max = 40, which is Sufficient to allow 
a bump extending over several windings while the radius 
return to the initial value at the end points, that is when 



-40 -30 -20 -10 0 10 20 30 40 
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FIG. 2: Mass and corresponding derivatives as function of an¬ 
gle, 0, for the parametrization in Ea. (l26ll . Top: m((^) (red), 
10 • m!{(f)) (green), 10 • m''{(f>) (blue), all in units of moR^. 

f.2 

Bottom contains potentials in units of that enter differ¬ 
ent quantization prescriptions: — (m" — /m^ (magenta), 

(m" —2^^)/m^(cyan), —K^/2-|-0.1(black), see Eas. lfT^ . (ITOl) 
and (flTl) . 

fx{4>) = fy{4>) ~ R. In numbers, 4>o = 27r means very 
fast variation over one winding. A given multiplum, £, of 
27r then implies the gaussian variation over 2^ windings. 
The geometric structure is illustrated schematically in 
fig© 

The parametrization in Eg. (1261) leads immediately 
from Eq.© to an effective mass given by 

m{4>) = moR^ 1^1 + . (27) 

The first and second order derivatives of the mass are 
then easily written down along with the potentials enter¬ 
ing the expressions for the different quantizations, that 
is Eqs. (O, (HU) and (HZl). We show their angular depen¬ 
dence in fig. HI 

The mass is constant for angles far away from the bulge 
on the circular helix. It exhibits a gaussian peak at the 
center with a width of 2(j)Q in accordance with Ea. (H71) . 
This is then the origin of a non-constant mass on the 
quantization, or equivalently the effect from quantal mo¬ 
tion in a one-dimensional curved space. 

The size of the effects is reflected in the variation of the 
derivatives shown in the same figure. Both first and sec¬ 
ond derivatives are rather small implying that the differ¬ 
ent quantizations at least qualitatively should reveal the 
same most important features. This is in spite of the ini¬ 
tial periodic parameterization where all such traces have 
disappeared since fx = fy and the square of the sine and 
cosine functions are added with equal amplitude. 
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FIG. 3: The absolute energies of the four lowest states of one 
particle on a wire parametrized in Eq. (1261) with a = 1 and 
</>o = drr. The energies are in units of f?/{R^mo). The upper 
panel shows the energies of Hemi and Hem^ where as the 
lower one shows the energies of /fgeo 


The differences in quantization are quantified by the 
three potentials at the bottom of fig. [21 For Hemi we 
have the potential appearing in Eg. (1131) after removal of 
the first order operator derivative. This term has an at¬ 
tractive minimum at the bulge and symmetric repulsive 
maxima in the tail of the bulge. The difference between 
HeMi and Hem 2 is comparable in size and almost flat 
over the bulge with two symmetric extremum points in 
the tail. 

The transverse-mode adiabatic approximation quanti¬ 
zation from Eq. dini) has two additional potentials com¬ 
pared to Hemi , that is one very similar to Eq. (fTO)) and 
the qualitatively different geometric potential propor¬ 
tional to the square of the curvature. The latter is al¬ 
ways attractive and of much larger magnitude, but oth¬ 
erwise with the same behavior as the mass itself, that is 
strongest in the center and vanishing outside the bulge, 
see fig. m 


B. Spectra 

We now solve the schrodinger equation numerically for 
the given parametrization for all the different choices of 
quantization described in section m We maintain the 
boundary conditions corresponding to a curve of finite 
length and fixed end-points, that '^{4>min maa:) — 

0. We discretize the angular space by choosing a grid, 
computing the finite element representation of the oper¬ 
ators, and diagonalizing in the corresponding basis. We 
then compare eigenvalues and eigenfunctions from the 
different quantization prescriptions. 

Inspection of the energies in tables U and m and fig¬ 


ure [3] reveal a clear pattern as all hamiltonians exhibit 
rather similar spectra, except for the ones with geomet¬ 
ric potential terms. We emphasize that the tables give 
absolute values for the ground state energies, and for bet¬ 
ter comparison the excitation energies are given for the 
excited states. 

Let us first consider the hamiltonians without the addi¬ 
tional geometric potential. Both absolute values and ex¬ 
citation energies for Hemi and Hem^ deviate from each 
other by at most 50% and the JWKB approximation al¬ 
ways assume intermediate values. All three energy sets 
approach each other with increasing excitation energy in 
accordance with an approach towards the validity of clas¬ 
sical physics. 

The simplest results are from the analytic JWKB ap¬ 
proximation which precisely is the n? spectrum from a 
one-dimensional square well. We know from Eq. (1241) that 
an overall energy scale is the square of the inverse length 
of the wire. This is seen as a constant ratio of 1.160 for 
JWKB energies of the same state. Applying this scaling 
from the JWKB approximation on the Hemi and He M 2 
spectra show the expected decreasing deviations for in¬ 
creasing excitation energies. These similarities already 
at the lowest energies are due to the relatively small co¬ 
ordinate variation of the effective mass. 

Including the geometric potential leads to substantially 
different spectra. This potential is overall attractive with 
constant curvature outside the bulge region in the cen¬ 
ter. The curvature and the attraction is smaller in the 
central region. This has strong implications on the re¬ 
sulting spectra where the lowest four states considered 
here all turn out to have negative energies. In figure [3] 
one sees how for larger a the states turns into two sets of 
degenerate states. This is because of the barrier in the 
geometric potential that shows up for larger a. 

Both JWKB and full solutions with the geometric po¬ 
tential have doubly degenerate ground states. This is 
due to the separation in two regions through the less at¬ 
tractive central peak of the geometric potential. There is 
room for bound states in each region, and no distinction 
between odd and even parity states. This is highlighted 
by the classically forbidden central region which is crucial 
for the simple JWKB solution. A better approximation 
would allow some tunneling into this barrier region and 
the degeneracies would be lifted as indicated by the en¬ 
ergies from the full solutions. 

The influence of the geometric potential is much larger 
than the variation between the Hemi and Hem^ spec¬ 
tra. This applies for both absolute and relative energies 
and independent of the bulge parameter a. In the limit 
of vanishing a the curvature approach a coordinate inde¬ 
pendent constant and the geometric potential becomes 
consequently constant as well. This amounts to a shift 
of all energies without any structural changes. The exci¬ 
tation spectra would thus approach each other. 
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TABLE I: The lowest four eigenvalues for different hamiltonians describing one particle on a wire parametrized in Eq. (I26II with 
a — 1 and 0o = dvr. The ground state energies are in absolute values in units of /{B?mo). The excited states are reported 
as excitation energies related to the corresponding ground state energy. . 


Hamiltonian 

Hemi 

Hem2 

JWKB 

Hgeo 

JWKB (including Vgeo) 

Ground state 

3.58 • 10“'‘ 

1.63 • 10“'* 

2.88 • 10"'^ 

-2.89 ■ 10"^ 

-2.76 ■ 10"^ 

1st excited 

7.18 • lO”"* 

1.02 • 10“® 

8.64 ■ lO"'^ 

8.62 ■ 10"® 

0 

2nd excited 

2.24 • 10“® 

2.42 • 10“® 

2.30 ■ 10"® 

5.07 • 10"® 

7.10 ■ 10"® 

3rd excited 

4.29 • 10“® 

4.37- 10“® 

4.32 • 10"® 

5.34 • 10"® 

7.10 • 10"® 


TABLE H: The same 

as in table U for a 

= 0.5 and 0o = dvr. . 


Hamiltonian 

Hemi 

Hem2 

JWKB 

Hgeo 

JWKB (including Vgeo) 

Ground state 

3.72 • 10“'‘ 

2.73 • 10“'‘ 

3.34 • 10""* 

-2.96 • 10"^ 

-2.86 • 10"^ 

1st excited 

9.01 • lO""* 

1.10 • 10“® 

1.00 ■ 10"® 

2.07 ■ 10"^ 

0 

2nd excited 

2.65 • 10“® 

2.70 • 10“® 

2.67 • 10"® 

3.14 • 10"® 

1.70 • 10"® 

3rd excited 

4.97 • 10“® 

5.03 • 10“® 

5.01 • 10"® 

4.76 • 10"® 

3.90 • 10"® 

n QH 


almost indistinguishable. 

The probability distributions 
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FIG. 4: Top: Ground state wave functions for the different 
choices of Hamiltonian, Eg. lfTSll (red), Eg. lfTTI) (green) Eo. (l25ll 
(black), Eo. (l20ll (blue). The constants in the parametrization 
are chosen to be a = 1 and 4>o = dvr. Bottom: The effective 
potentials for He M2 (red) and Hgeo (blue) for the curve in 
Eo. (l26ll with a = 1 and (j)o = dvr 


C. Eigenfunctions 

The complete picture requires spectra supplemented 
by corresponding eigenfunctions. The ground state wave- 
functions for the parameter choice a = 1 and = drr are 
shown in fig. 0] for the different hamiltonians. All the so¬ 
lutions must vanish at the end of the wire, since this is 
the imposed boundary conditions. 

Let us again first consider the cases without geomet¬ 
ric potential. The ground state solutions are all sym¬ 
metric around the center where the radius is largest and 
the largest probabilities occur as peaks. The solutions to 
HEM^ in Eq. (HU and the JWKB result from Eq (1301) are 


are rather broad and extending far beyond the bulge in 
the central region of the curve. The solution for Hem 2 in 
Eq. (1151) is rather similar although distinguishable with 
a central peak, slightly higher and correspondingly nar¬ 
rower. 

The difference between the Hemi and He M 2 solutions 
arises from the potential in Eq. (fTOl) . ; which 

is plotted in fig. [H This potential is zero at both ends, 
becomes repulsive when moving towards the center from 
either side of the potential, and finally it turns attractive 
in the center. The wavefunction exhibits a sharp upwards 
turn when the attraction is felt with the result of a larger 
maximum than for Hemi ■ 

Inclusion of the geometric potential changes qualita¬ 
tively all the derived solutions. The full numerical solu¬ 
tion to Hgeo is still symmetric but with an almost vanish¬ 
ing minimum in the center and two prominent maxima 
on either side. This behavior is a direct reflection of the 
properties of the geometric potential shown in fig. 0] and 
contained in Hgeo- The JWKB solution now has two sep¬ 
arated classically allowed regions, of course provided the 
energies are below the barrier in the center. The simplest 
JWKB solutions are then zero in the central forbidden 
region as seen in fig. |4l The deviation from the full so¬ 
lution is therefore very striking but uninteresting since 
the probabilities at the same time are very small. An 
improved JWKB solution could be designed by use of an 
exponentially decreasing wave function to describe tun¬ 
neling into the barrier. 

The first excited states are shown in fig.[5]for the differ¬ 
ent hamiltonians. The boundary conditions of zero at the 
end points of the wire are maintained in all cases. Now 
a node appears in the center, and all wave functions are 
of odd parity. The solutions without geometric poten¬ 
tial are almost indistinguishable for the Hemi, Hem 2 ^ 























2 ; 



FIG. 5: First excited states for the different choices of Hamil¬ 
tonian, Eg. lfT^ (red), Eg. lfTTI) (green), Ea. (l^ (black). The 
constants in the parametrization are chosen to be a = 1 and 
00 = drr. 

solutions and the related JWKB approximation. They 
all resemble the sine wave functions for the first excited 
state of a particle in a square well potential. Now the 
central attraction for He M 2 is much less effective due to 
the required node for 0 = 0. 

Again including the geometric potential changes the 
solutions substantially, although much less than for the 
ground state. The odd parity characteristics with a cen¬ 
tral node is maintained but now with a small oscillatory 
modulation by tunneling into the barrier. The two re¬ 
gions of large probability are pushed further away from 
the center by the potential than the Hemi and He M 2 
solutions. The slope of the wave functions is also not as 
steep across the central region. 

The corresponding simplest JWKB solution is trying 
to mimic this behavior in the classically allowed regions. 
The central forbidden region has a constant probability 
of zero. The only difference between first excited and 
ground state is that the ground state is even, and the first 
excited state is odd. The absolute values corresponding 
to the probability distributions would be identical. 

The dependence of the eigenfunctions on the parame¬ 
ters of the gaussian central bump is intuitively clear with 
the detailed knowledge we accumulated from the investi¬ 
gated set. The effects from varying the two parameters, 
a and 0Oi seems to be very different. The value of aR 
is directly the radial extension of the bump beyond ini¬ 
tial helix radius. It is therefore clear that increasing a 
from zero must increase the change of the solutions from 
the a = 0 solutions where both mass and curvature are 
constants. However, the qualitative behavior of even and 
odd parity is maintained, and with the related maxima 
or nodes at the center. 

The 00 variation appears to be very different but the 
effects are actually rather similar. Large values imply 
a slow variation of the radius of the helix, and as such 
small influence on the wave functions beyond a possible 


FIG. 6: Top: Ground states for the different choices of Hamil¬ 
tonian, Eo. IT^ (red), Eg. JTTI) (green), Eo. ll^ (black). The 
constants in the parametrization are chosen to be a = 0.5 and 
00 = dvr. Bottom: The effective potentials for H'_BM 2 (red) and 
Hgeo(blue) for the curve in Eg. (1261) with a = 0.5 and 0o = Jtt 


scaling from a different average radius. Small 0o values 
imply fast variation over few windings. Both mass and 
curvature would then vary much faster as well, and the 
different quantizations would be very different. 

Thus, large a and small 0o lead to large variation in 
mass and curvature and consequently the different quan¬ 
tization prescriptions would deviate more and more from 
each other. This would be particularly prominent in com¬ 
parison with use of the geometric potential. It is worth to 
emphasize that it is not obvious which quantization pro¬ 
cedure is most correct for these one-dimensional cases. 

On one hand the simple JWKB approximation pro¬ 
vides a very accurate match with the Hemi solutions. 
However, this assumes that no potential is necessary to 
confine the particle to the one-dimensional wire. On the 
other hand, some geometric potential combined with an 
appropriate kinetic energy operator would directly de¬ 
liver the hamiltonian. Unfortunately, this assumes a 
computational scheme to obtain a reliable potential, and 
the lowest order curvature dependent potential is not 
accurate for helix like periodic structures with strongly 
varying effective mass. 


V. STRETCHED HELIX 

In contrast to the previous section we shall here inves¬ 
tigate asymmetric helix deformations. We design here 
two stretching parameterizations originating from very 
different assumptions. We first describe these parame¬ 
terizations, and in the following subsections we present 
results for spectra and eigenfunctions for the different 
quantization descriptions. 
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A. Parameterizations 

We want to study the non-trivial cases where both 
mass and curvature are monotonously varying with the 
coordinate. Instead of symmetry we choose an increasing 
stretching along the symmetry axis of the helix, that is 

/.(</>) = /yW ,/.(</>) (28) 

where a is a constant. This curve is circular in the x — y 
plane and the distance between the windings increase 
with the angle, </>, see fig. [71 The effective mass is simple, 
that is 


= moE? (l + , (29) 

which allow analytical integration of the square root in 
Eq. (1751) . and therefore a fully analytic JWKB-solution. 
Only first and second derivatives are finite and expansion 
in higher order derivatives are more likely to converge 
than for a periodic structure. 

Explicitly we get the bound state wave function given 
by Eq. (1251) with 


L{(j)) = R 




+a4) 



1 

4a^ 



where L{(f>max) is the length of the wire. 

The two hamiltonians Hemi and He M 2 differ by the 
second derivative of the inverse mass, see Eq. (fT6|) . We 
can find a parametrization where this difference van¬ 
ishes. The assumption of identical circles in the x — y 
plane, fx = fy = ^, gives an effective mass from Eq. o, 
m{(j)) = too-R^(1 -|- [fzY)- If we therefore assume that 
He Ml = He M 2 then 1/(1 + (/z)^) should be a first order 
polynomium i 0, or equivalently m{(j)) = where 

Co and Ci are constants. 

Integrating //(</) to find fz{4>)-: we get the parametriza¬ 
tion: 


/.(</) 

/.(</) 


fyW = 1 


arctan 

— 2co—2ci(ti+l 

2-y/-(oo+ci<^-l)(co+ci0) 


2ci 


2RiJ- (cq -b Cl/) - 1) (cq -b Cl/)) 


The curves in Eqs. (l29l) and (1311) are different types of 
monotonous deformations in the z-direction. The result¬ 
ing deformed helices are both shown in fig. [71 

The first and second order derivatives of the mass are 
then easily calculated. We show their angular depen¬ 
dence in figs. ISI and [5] and as well the expressions in 
Eqs. (O, (USD and (II3 that enters the expressions for 


\ 

\ 

; 

/ 

i 



FIG. 7: Schematic drawings of two deformed helices with 
stretched and squeezed parametrizations Eqs. dSSl) and m- 
The stretching and squeezing parameters are chosen to be 
a — 0.1, and cq = ci = 0.01. 



FIG. 8: Mass and corresponding derivatives as function of 
angle, 4>y for the parametrization in eg. (1281) . Top: (red), 

■m!{(j)) (green), ■m”{(j)) (blue), all in units of moR^. Bot- 

f.2 

tom contains potentials in units of that enter differ¬ 

ent quantization prescriptions: — (m" — §^^)/m^(magenta), 
(m” — 2^22_)/m^(cyan), —k^/ 2(black), see Eos. (1131) . (1161) and 

Ha- 


tire different quantizations. The two parameterizations 
can be viewed as stretching and squeezing, respectively. 
The structure variations in these quantities are therefore 
rather similar, except that they appear at small or large 
/>, respectively as seen in figs. |5| and IH 

The mass itself increases quadratically or decreases in¬ 
versely proportional with (f>. This smooth dependence is 
then the origin of the effects of this type of non-constant 
mass on the quantization. The different combinations 






























10 



0 5 10 15 20 25 30 35 40 


FIG. 9: Mass and corresponding derivatives as function of 
angle, 4>, for the parametrization in eg. (1311) . Top: m((/>) (red), 
(green), -m"{(f)) (blue), all in units of moR^. Bot- 
tom contains potentials in units of that enter differ¬ 

ent quantization prescriptions: — (m" — |(magenta), 

(m" — 2-^^)/m^(cyan), — j^K^/2(black), see Eas. dTSll . (ITOl) 
and (fTTl) . 

are then very trivial as for example both m" and m''^ jm 
are constants. The difference between Hemi and Hem 2 
therefore quickly vanishes with (/) in the first case while 
by construction identically equal to zero in the last case. 
The different combinations exhibit the opposite behavior 
with (j). The additional potential in Eq. (l20l) is very small 
for both parameterizations. On the other hand, the ge¬ 
ometric potential is again very decisive with prominent 
minima at either small or large values of (f). 


B. Spectra 

In tables IIIII and IIVI we show the lowest four states of 
the excitation spectra of each hamiltonian. The Hemi 
and He M 2 energies for both parameterizations are very 
close to the corresponding JWKB spectra where the dif¬ 
ferences almost disappear for the highest excited states. 
The deviation is largest between the absolute values of 
the ground state energies of the Hemi and HEM 2 hamil- 
tonians. The largest differences between the parameter¬ 
izations can be removed by the length scaling which is 
cleanly expressed by the analytic JWKB expression in 
Eq. 

Inclusion of the geometric potential changes the spec¬ 
tra substantially. This potential is attractive and able to 
support one or two bound state with negative energy, re¬ 
spectively for the two parameterizations. These features 
are also found in the corresponding JWKB spectra. Fur¬ 
thermore, the JWKB excitation energies are approached 



FIG. 10: Ground states for the different choices of hamil¬ 
tonian, Eg. lfTSl) (red), Eq.f lfTTIl (green), Ea. (l25l) (black) 
for the parametrization in Eg. (12811 . The constants in the 
parametrization are chosen to be a = 0.1. 


for higher excitations. The length scaling is not appro¬ 
priate here, since the wave functions are pulled into the 
attraction regions which eliminate the importance of the 
finite size confinement. 


C. Eigenfunctions of stretched helix 

The stretched helix is parametrized by Eq. (E51) where 
we choose again the stretching parameter a to have the 
value a = 0.1. The parametric angle interval for the curve 
is (p varying between (pmin = 0 and (pmax = 20. The 
ground state wavefunctions for all the different hamilto- 
nians are shown in fig. 1101 The solution to Hemi and 
He M 2 and the corresponding JWKB result all exhibit a 
single maximum shifted from the center at = 10 to¬ 
wards higher values of (j). These wave functions are very 
similar, although that of Hemi has the maximum shifted 
a little more than the almost indistinguishable results for 
He M 2 and the JWKB approximation. These shifts are 
towards higher values of the effective mass and all due to 
the corresponding increase with ip. 

The geometric potential has a very strong effect. The 
ground state wavefunctions still only have one peak, but 
now shifted towards smaller values of p, where the curva¬ 
ture is larger and the attraction therefore stronger. The 
corresponding JWKB solution is similar with one peak at 
roughly the same position as the full solution. However, 
the classically forbidden region of VgeoiP) > E starts at 
p = 10, after which the JWKB wavefunction is zero by 
definition. 

The wavefunctions for the first excited state are shown 
in fig. 111! where the necessary node is the prominent fea¬ 
ture. Again almost quantitative agreement within the 
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TABLE III: The lowest four eigenvalues for different hamiltonians describing one particle on a wire parametrized in Eq. (1281) 
with a = 0.1 and 0 varying between 4)min = 0 and (pmax = 20. The ground state energies are in absolute values in units of 
/{B?mo). The excited states are reported as excitation energies related to the corresponding ground state energy. 


Hamiltonian 

Hemi 

Hem2 

JWKB 

Hgeo 

JWKB (including Vgeo) 

Ground state 

2.72 • 10“® 

1.57- 10“® 

2.29 ■ 10"® 

-1.51 ■ 10"^ 

-6.01 ■ 10"® 

1st excited 

7.00 • 10"® 

6.59 • 10"^ 

6.85 • 10"® 

1.69 ■ 10"® 

7.46 • 10"® 

2nd excited 

1.84 • 10“^ 

7.72 • 10“^ 

1.84 • 10"^ 

2.63 • 10“® 

1.70 ■ 10"^ 

3rd excited 

3.43 • 10”^ 

9.31 • 10”^ 

3.44 • 10"^ 

4.11 • 10"® 

3.19 • 10"^ 


TABLE IV: The same as table Hill for the parametrization Eq. m with Co = Cl = 0.01 and 4> varying between 4>min = 0 and 
(pmax — 40. 


Hamiltonian 

Hemi 

Hem2 

JWKB 

Hgeo 

JWKB (including Vgeo) 

Ground state 

3.93 • 10-'‘ 

3.93- 10-'‘ 

4.22 • 10"^ 

-8.99 ■ 10"® 

-6.99 ■ 10"® 

1st excited 

1.26 • 10"® 

1.26 • 10“® 

1.27- 10"® 

6.23 ■ 10"® 

5.25 • 10"® 

2nd excited 

3.37- 10"® 

3.37- 10"® 

3.38 • 10"® 

9.23 • 10"® 

7.38 • 10"® 

3rd excited 

6.32 • 10“® 

6.32 ■ 10“® 

6.34 • 10"® 

1.17 • 10"® 

9.77 ■ 10"® 



FIG. 11: Excited states for the different choices of hamil- 
tonian, Eg. JTSl) (red), Eo.l' lfTTIl (green), Eg. (1251) (black) 
for the parametrization in Eo. (l28l) . The constants in the 
parametrization are chosen to be a = 0.1. 


two groups of results, that is between the Hemi, Hemz 
and related JWKB results, and between the Hgeo and 
corresponding JWKB results. The added geometric po¬ 
tential changes the quantitative behavior by moving the 
peaks towards smaller (()-values where the attraction is 
largest. 


D. Eigenfunctions of squeezed helix 

The squeezed helix, which is parameterized by 
Eq. (1511) . is designed to give Hemi = Hem^- The ef- 



FIG. 12: Ground states for the different choices of hamil- 
tonian, Eo. JTSl) (red), Eo.(' (fTT1) (green). Eg. (1251) (black) 
for the parametrization in Eo. (l28|l . The constants in the 
parametrization are chosen to be co = 0.01 and ci = 0.01 


fective mass decreases with </> as m((/)) = where 

we choose cq = ci = 0.01 and a curve parametrized by 
(f> varying between (pmin = 0 and ipmax = 40. We show 
the ground state wavefunctions in fig. [121 They are al¬ 
most a left-right reflection of the stretched wave func¬ 
tions in fig.dni The overlapping wavefunctions of Hemi 
and Hem 2 rise quickly from zero to a single maximum 
at (() = 9, before they linearly fall off to 0 at (pmax- The 
JWKB solution is similar and rises from 0 to a maximum 
at (j) = 12, before it falls off to zero at (j>max- 

Again the geometric potential moves the peak to the 



























12 



<P 


FIG. 13: First excited states for the different choices of 
hamiltonian, Eg. dlSIl (red), Eg.f dllll (green), Eo. d25l) (black) 
for the parametrization in Eo. d28ll . The constants in the 
parametrization are chosen to be co = 0.01 and ci = 0.01 


opposite end of the allowed 4> interval, that is to about 
(/) = 32. The corresponding JWKB solution is similar 
with a maximum at roughly the same </> value. The de¬ 
crease is steeper towards the classically forbidden region 
for (/) < 22. Thus the picture is that the geometric po¬ 
tential move the solutions to the large curvature region, 
which is the opposite of the large effective mass region 
where the Hemi and Hem 2 solutions are peaked. 

The first excited states of the same configuration are 
shown in figure 1131 They all now have the required node 
for an excited state. The overlapping solutions from 
He Ml and He M 2 rise to their first maximum at ^ = 4, 
then they have a node at (j) = 11, and then a smaller 
minimum aX (p = 23. The related JWKB solutions is 
similar with slightly shifted extremum points. The ge¬ 
ometric potential leads to a wavefunction with a broad 
peak at the center, and a smaller minimum at the posi¬ 
tion where the ground state wavefunction is peaked. This 
is required by orthogonality. The corresponding JWKB 
solution is similar but has as usual to vanish within the 
classical forbidden region for </> < 10. 


VI. DISCUSSION AND OUTLOOK 

We start with two different approaches to a system of 
a single particle trapped in an effective one dimensional 
trap. The first approach is to build a classical descrip¬ 
tion of the system, and through that find an appropriate 
quantization. Because we allowed this one dimensional 
trap to have a changing curvature, this quantization step 
is not trivial and we show two equally valid choices, which 
differ only by a potential term. 

The other approach starts from a quantum mechan¬ 
ical description in three dimensions, and then through 
a transverse-mode adiabatic approximation reduce to an 
effective one dimensional model but now with an extra 
so-called geometric potential. This potential is attractive 
and given as proportional to the square of the curvature. 
We then investigate three different perturbations of a he¬ 
lix, and calculate the wavefunctions and energies of the 
different hamiltonians. 

Monotonous deformation results in monotonous effec¬ 
tive mass and curvature dependence on the coordinate. 
However, these two key quantities behave differently 
and lead to opposite effects on the quantized solutions. 
Specifically, increasing curvature leads to increasing at¬ 
traction along the wire, and therefore the ground state 
wavefunctions would peak at this end. Increasing effec¬ 
tive mass also tend to move the largest probability in the 
same direction of large mass. This implies that the differ¬ 
ent quantization prescriptions in this case of monotonous 
helix deformation produce very different results. 

The more periodic type of helix deformation leads to 
more similar quantized results although still with sub¬ 
stantial differences. The periodic nature of a helix pro¬ 
hibits that the confining potential is obtained by a con¬ 
verged Taylor expansion in terms of the parametrizing 
one-dimensional path. It also strongly indicates the same 
problem with a quantization obtained by forced, explicit 
symmetrization of a non-hermitian hamiltonian. The 
only hamiltonian without these problems has the inverse 
effective mass between the two derivatives in the kinetic 
energy operator. 

The authors acknowledge inspiring conversations with 
J. Stockhofe and A. Rauschenbeutel. This work was sup¬ 
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